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Abstract 

We consider a symmetrical star-shaped network, in which bandwidth is 
shared among the active connections according to the "min" policy. Start- 
ing from a chaos propagation hypothesis, valid when the system is large 
enough, one can write equilibrium equations for an arbitrary link of the 
network. This paper describes an approach based on functional analysis 
of nonlinear integral operators, which allows to characterize quantitatively 
the behaviour of the network under heavy load conditions. 

1 Model description 

Consider a network comprising N links, where several data sources establish 
connections along routes going through these links. The main concern here 
is about the policies that can be used to share the bandwidth of the links 
between active connections, and the effect of these policies on the dynamics of 
the network. 

In this paper, the network is star-shaped (see Figure [lj and all routes are of 
length 2, which is a reasonable model for a router. Each star branch contains 
two links ("in" and "out") and each route is isomorphic to a pair of links. 

Let r = denote a route on links i and j, and 1Z be the set of all 

possible routes (with cardinal N 2 /2). Connections are created on r according 
to a Poisson process with intensity 2X/N, so that the total arrival intensity on 
each link is A. A connection lasts until it has transmitted over the network its 
data, the volume of which follows an exponential law with mean v. Each link 
i, 1 < i < N, has a bandwidth equal to 1 and its load is p = Av. 

The state of the system at time t £ IR is given by the number of active 
connections on each route (c ( r N) (t), r e TZ) = X (N \t). The vector X (N) {t) is in 
general Markovian and 
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Figure 1: The star-shaped network 



is the total number of active connections on link i, VI < i < N . 

It is now necessary to describe how bandwidth allocation is achieved. The 
"max-min fairness" policy, popular in telecommunication models, being too diffi- 
cult to be studied rigorously, this article focuses on the "min" policy, proposed by 
L. Massoullie and J. Roberts [5], in which a connection on gets bandwidth 

min (wWYYfW))- (L1) 

This allocation clearly satisfies the capacity constraints of the system, and 
can be shown to be sub-optimal with respect to max-min fairness. Its invari- 
ant measure equations are however too complicated to be solved explicitly. An 
efficient method in such situations is to study the model in so-called thermody- 
namical limit, using mean field analysis. 

In order to study the stationary behaviour of the network as N — > oo, the 
idea is to assume the conditions of chaos propagation, under which any finite 
number of links tend to become mutually independent. In this paper, this 
hypothesis will be considered as a heuristic, to be proved in further studies. 
Some rigorous studies of this type can be found in Vvedenskaya et al. [5] and 
Delcoigne et Fayolle [SJ. 

From now on, p < 1 and the system is assumed to be in stationary state 
X = (c r ,r 6 TV). For all k > 0, the empirical measure of the number of links 
with k connections is 

(N) dcf 1 \ ^ fl 

l<i<N 
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Symmetry considerations imply that, for all i < N, P(X 4 ' JV) — k) = ~E,a k N) . 
Besides, from the chaos propagation hypothesis, a law of large numbers is as- 
sumed to hold for a k N) : 

a fc = lim ajW= lim P(A 4 (iv) = k) = W{X = k). 

N—t-oo N—>oo 

The a k s, traditionally named mean field, drives the dynamics of the system. 
The following notation will also prove useful: 

aC'^Elf = ^ fcojW , a = ^ fca fe . 

fc>0 fe>0 

A heuristic computation (detailed in [3]), yields the following equations: 

a k+ iu k+1 = paa kl Vfe > 0. (1.2) 

While (|1.2p resembles a "birth and death process" equation, it is in fact highly 
non-linear, due to the form of u k and of a. 

The purpose of this paper is to show how the asymptotic behaviour of the 
system (as p — > 1) can be derived from the analytical study of the generating 
functions built from (|1.2[) . This work is a part of the wider study [3], which also 
gives ergodicity conditions for any topology under the min and max-min policies, 
shows how equations like (|1.2[) are derived (also in the case where routes are 
longer than 2) and presents comprehensive numerical results. 

The main byproduct of Theorem l3.3l is the following asymptotic expansions, 
valid in a neighborhood of p = 1. 



lim p ak ~ (1 — p)B exp 

k—¥oo 



(1 - p)A 



where A and B are non-negative constants. Moreover, if c(z, 1) and v(z, 1) are 
the solutions of the system of differential equations ()3.5[) , then A can be written 
as follows: 

poo 

A = I c(z,l)dz= lim zv'(z, 1) » 1.30. 



Since this system is numerically highly unstable, it has proven difficult (with 
the "Livermore stiff ODE" solver from MAPLE) to get a better estimate for A. 



2 An integral equation for the generating func- 
tion 

Let C(r) (resp. T>(r)), be the circle (resp. the open disk) of radius r in the 
complex plane. 
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Let a : z — » a(z) be the generating function, a priori defined for z in the 
closed unit disk 

/ \ dof k 

a(z) = 2_^a k z . 

k>0 

Denoting a' the derivative of a, (|1.2p can be rewritten as 

a k+1 u k+ i = pa'(l)a k , Vfc > 0. (2.1) 

Lemma 2.1. 

(a) If has a probabilistic solution, then, necessarily p < 1, and 

lim a k p~ k = K(p), 

k—>oo 

where K(p) is a positive constant, bounded^ p < 1. 

(b) The function a satisfies the nonlinear integral equation 

a' (1)(1 - pz)a(z) = -J- / a(iv)a(-) - — , (2.2) 

2i7r JC(r) \U)/[l-U)y 

where \z\ < p^ 1 et r is an arbitrary positive number, with 1 < r < p^ 1 . 



Proof. When (|1.2p has a probabilistic solution, necessarily a = lim„ i _ i . 0O u m . 
Therefore, Ve > 0, there exists a number M(e) > 0, such that 

p <^l+l<^ i Vm>M(e), 
a m 1 - e 

which implies that p < 1. Moreover, under the same existence hypothesis, one 
can write 

k 

a k p- k = a T[ _ a , (2.3) 

where = S m >£( m ~ 

When £ — > oo, the convergence of the product in (|2.3|) is equivalent to the 
convergence of the series of general term Dg, which holds since 

D t < a ( m ~ Q{P + e )™ W > M ( £ )' 
so that ^ f Di behaves like 

<>0n>0 v ^ ' 

Point (a) of the lemma is proven. Point (b) is an application of (a) and of 
an integral representation used by Hadamard, recalled below (see e.g. [7]). 
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Let a and b be functions 

a(z) = 2J a n z n , b(z) = 2J &n2™, 

n>0 n>0 

analytic in the respective disks T>{R) and T>{R'). Then the function 

c(z) = y^a„6 K z" 

n 

has a radius of convergence greater than RR' and has the integral form 

1 f , v . / z\ du> 



c(z) = — / a (w)&(-)— , 

Z17T //■ \W/ CO 



c 

where C is a closed contour containing the origin, and on which \co\ < 
R, - <R'. 

From (a), a has a radius of convergence at least equal to p^ 1 : this property, 
used in Hadamard's formula, leads directly to (|2.2[) . □ 



The following Lemma provides a finer description of a(z). 
Lemma 2.2. The function a is meromorphic and can be written as 

oo 

«(*) = £^r— • (2-4) 

t-^ p 1 — z 

1=1 ' 



Moreover, for any sequence of circles C(R n ), such that 
(1 + e)p-" <R n <(l- e)p-<" +1 \ < e < 



1 + p- 

one has \a(R n )\ = o(l) as n — > oo. 



Proof. The integral equation (|2.2p allows for the analytic continuation of a in 
the whole complex plane. Indeed, from point (a) of Lemma l2.11 a is holomorphic 
in V{p~ 1 ) and its first singularity is a simple pole at the point z — p^ 1 . An 
application of Cauchy's theorem to the integral in (|2.2[) leads to 

7T- / =J-f + Residue^- 1 ), 



2in Jc(Ri) 2ilT JC(r) 

where Ri is defined above. Since a(uj) and a(z/u>) are analytic in the regions 



< P -\ 



the integral is a function of z analytic in the ring-shaped area p^ 1 < \z\ < p~ 2 . 
The same holds for a, thanks to the left-hand side of (|2.2[) . Equation (|2.4p 
follows by recurrence. 
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The second part of the lemma is obtained by application of Cauchy's theorem 
to the integral in (|2.2[) along the circle C(R n ): 



a'(l)(l - pz)a(z) - £ jfi*^ = a(«) a(£) ^ 



da; 



One can make the analytic continuation of the above equality, the left-hand 
side of which is analytic in \z\ < R n , letting z reach the circle C(R n ) along a 
simple curve which avoids the polar singularities = l,...,n. Bounding 

the modulus of the integral yields the inequality 



( P R n - l)M(Rn) < A n sup (M(Ri)) + DM{R n ) 

2<n— 1 



\duj\ 



c(i?.„) I 1 ~ w l 2 ' 



where D is a positive constant and A n is bounded Vn < oo. By induction, 
limn^oo M(R n ) — follows easily, and the proof is completed. □ 

Combining (|2.4p with ()2.2() , a straightforward computation of residues gives 

OO OO _j 

a'(l)(l pz) Y —2— = Y - ° iajf !, ^. (2.5) 

L -tp- % -Z . Z -^(p- ( -' l +^-z)(l-p-3) 2 V ' 

Let a : z — > a{z) be the generating function 

a(z) = a k+1 z , 
fc>o 

defined for z in a bounded domain of the complex plane, including the origin. 
Using point (a) of Lemma 12.11 and (|2.4p , the following relations hold 

ai = pK(p) 

(2.6) 

a k - X = p k a(p k ), Vfc > 1, 

Since a < 1 — p, the function a is thus analytic in the disk T>(p). Identifying 
the coefficients of the power series in z in (|2.5p . one gets 

It follows easily by recurrence that the a^s are of alternate signs, with a\ > 0. 
Let 
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Hadamard's formula, when applied in (|2.7[) , implies the integro-differential 
equation 

a'(l)[a(pz)-a(z)] = ^- [ a(w)/(-W (2.9) 
2^7^ J c(r) \ljJ 

valid in the domain {r < p, \z\ < p~ 2 }- 

Taking the second form for / in (|2.8p , which in fact converges in the domain 
{p < 1, $t(z) < 0}, and applying Cauchy's theorem to the integral in (|2.9|) . one 
obtains the functional equation 

a'(l) [a(pz) - a(zj] = a(pz)b(z), \z\ < p, p < 1, (2.10) 

where 

oo 

6(z) d = f z^jy +1 a(^ +1 z). 
i=i 

From (|2.10p . let's make now the analytic continuation of a in the nested 
disks U(p~ n ),n > 1. It appears that a has no singularity at finite distance, and 
consequently is an integral function. From the general theory [7!, it follows that 
a is completely characterized by its zeros and its order at infinity. 

Let zq be an arbitrary zero of a. From (|2.10p again, a(p~ l zo) = and the 
zeros of a form families of points in geometric progression with parameter p^ 1 . 
It suffices to determine the zeros of smallest modulus, but, alas they do not 
have any explicit form and numerical schemes are highly unstable. However, 
from ([231), K 1 ) = "'(!)> so that (i2~T0|) implies a(l) = together with 

a(p-' l )=0, \fi>0. (2.11) 



3 On the asymptotic behaviour around p = 1 



In order to assess the practical value of the "min" policy, it is important to evalu- 
ate the system behaviour in heavy traffic conditions. The numerical calculations 
in [3] show that the distribution of any Xi is modal, which is not common in 
known models. 

In this section, it will be convenient to consider p not only as a parameter 
but as a plain variable. Therefore, in all quantities of interest, p will appear as 
an explicit variable, e.g. a(z,p), f(z,p) or ak{p)- 

The fundamental ideas of the analysis will be given after the next lemma, 
which proposes a scaling — likely to be the only interesting one — for the function 
a(z,p). 



Let 



Ck(p) ' 
c(z,p) 



-ai(p) 



(logp) 3 a'(l, j o)' 

ef ~ a k(p) 

= (logp)3a'(l,p)e(p) fe! 

oo 

dcf i \ k 

= 2_^c k+ i{p)z . 



Vfc > 1, 



(3.1) 
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The reader will easily convince himself that the factor (1 — p) 3 arises rather 
naturally; however the factor — log 3 p has been chosen here, since it provides 
more compact formulas in the forthcoming results. 



Lemma 3.1. Let 



v(z,p) 



(logp) ; 

2l7T 



C(l) 



c(u,P)f[-,p)du, \ Z \<P 



The following functional relations hold: 

'a(z,p) = a 1 (p) c(z£(p),p), 
c(0,p) = l, 

c(z, p) = c(pz, p) [1 + log pv(z, p)] , 



(3.2) 



c(pz, p) [1 + z(log pf ip l+1 c{p l+1 z, p)] . 



The coefficients Cfe(p), k > 1, are of alternate signs and the function c(z,p) 
has the following properties. 

(a) there exists only one solution c(z, p) of H3.2\) , which is integral with respect 
to z and bi-analytic in (z, p) in the region < p < 1. 

(b) Define 



3=1 3 J ° 



log(l — u)du 



V|t| < 1, 



denoted by some authors as dilog(l — t). Then c{z, 1) exists and satisfies the 
integro- differential equation 



dc(z,l) _ c(z,l) 
dz 2iw 

which rewrites in the form 

dc(z, 1) 



C(r 



c(u;,l)g(^jdu, \r\ = 1, V|z| < 1, 



c( Z ,i) r C ( W ,i)iog(i) 



dz 

which is equivalent to the non-linear differential system 

dc(z,l) 



diu, 



(3.3) 



(3.4) 



dz 



+ c(z, l)v(z, 1) = 0, 



with initial conditions 

w(0,l) = 0, 



d 2 v{z,l) dv(z,l) 
dz 2 dz 

dv(z, 1) 
dz \ z =o 



(3.5) 



c(z,l), 



1, c(0,l) = l. 



(c) Moreover, c(z, 1) is analytic in the open complex plane, except at a neg- 
ative real point q, and c(z, 1)^0, Vz^qUoo. 
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Proof. The first three equations in p. 21) follow directly from the definition (|2.7p 
and (|2.10[) of the coefficients a,k(p), the fourth one coming from the analytic 
continuation of (|2.9[) . Existence and uniqueness are simple consequences of the 
convolution equation (12.71) . 

The properties relative to the morphology of c(z, p), p < 1, are more intricate. 
First, the reader will notice that (|3.3p can be obtained rigorously from (|2.7|) 
or (|2.8p . but not from (|3.2p ! Then, there is a phase transition when p = 1. We 
will return to this topic in Section 01 □ 

It is interesting to note that the function w(y) = v(e y ) + 1 satisfies the 
so-called Blasius [1] third-order differential equation 

w"'(y) + w(y)w"(y)=0, 

which arises in hydrodynamics to describe the stationary evolution of a laminar 
boundary layer along a flat plate! The explicit solution of this equation is still 
unknown, albeit it has been studied by many authors over the last decades (see 
e.g. 0). 

Starting from Lemma l3.il it is now possible to sketch the main ideas of the 
proposed method. The Gordian knot amounts to the evaluation of £(p), defined 
in (|3.1[) . This can be done via the anchoring equation 

c(a P ),p)=0, (3.6) 

which follows from 12 .111 and from the first equation of p.2[) . From the structure 
of the third equation of (|3.2[) . it appears that the smallest positive solution of 
c(it, p) = satisfies 

— ^ = (logp) 2 £^+ 1 C (V+\p), 
ulogp ^ 

the right-hand side of which is an analytic function, bounded in any compact 
set Vp < 1 : as p — > 1 , necessarily u — > oo and all positive zeros of the anchoring 
equation are sent to infinity. The key is to find the asymptotic behaviour in z 
of the various functions, in the cone < z < U(p), which contains £(p): in this 
cone, c{z,p) is close to c(z, 1) — in some sense — and 

v(z,p)*w(z,p) d = f log 2 p^V+ 1 C (^+ 1 I l). (3.7) 

Since w(z, 1) has a logarithmic behaviour, £(p) can be obtained by direct 
inversion. The sketch of the proof is outlined below. 

We will need the Mellin transform (see e.g. [3]) of c(z, 1), defined as 

POO 

c*(s) = / x s ~ 1 c(x)dx. 
Jo 

The behaviour of c(z, 1) in the region 5ft(z) > 0, given in Section 21 implies 
the existence of c*(s, 1), Vs, Sft(s) > and of all the moments of c(z, 1) on the 
positive real axis. 
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Lemma 3.2. The junction w(z,p) defined by \3. 7)) admits, Vz, 5R(z) > 0, £/ie 
asymptotic expansion 



w(z,p) = c*(l, l)log(pz) - 



dc*(l,l) 
ds 



+ (log p) 2 [(log p) + p)] + 0(z- d ) 



(3.8) 



where d is an arbitrary positive number, $ and VP are fluctuating functions of 
small amplitude for p =/= 1, which vanish for p = 1. 



Proof. From the Mellin transform inversion formulas, 



tr+ioo 



2-V(« + l,l) 



logp 



1 - p s 



ds, We (-1,0). 



(3.9) 



Let s„ = 2m7rlog 1 p, V?i G Z. Cauchy's theorem can be applied to 
by integrating along the vertical line 3J(s) = d > 0, so that 



w 



n£Z 



c*(s„ + l,l)log(pz) 



dc*(s n + 1,1) 



+ (3.10) 



The series above is equal to the sum of the residues, taken on the vertical line 
5R(s) = 0, and is uniformly bounded, Vp < 1. Indeed, an integration by parts 
gives the inequality 



|c*(s n + 1,1)1 = f 
Jo 



y Sn c(y, l)dy 



< 



T{ Sr , 



1) 



r(s r , 



k + l) 



y k \c^{y,\)\dy, 



where is the fc-th derivative of c[z, 1) and T(x) is the usual Eulerian function. 
Since 

r(s„ + i) 



< 



log"' P 

k\ ' 



r(s„ + k + 1) 

the proof of (|3.8p and of the lemma is concluded. □ 
Theorem 3.3. For some real number d > I, the following expansions hold. 

1 <91ogc*(l,l)" 



p£(p) = exp 
a'(l,p) - 



logpc*(l,l) 
1 



ds 

0((logp) d ). 



;i + 0((logp) d )),(3.11) 
(3.12) 



(log p) 2 c* (1,1) 

Proof. From (|3.6p . £(p) is solution of the equation in x 

1 + log pv(x,p) = 0. 

A deep analysis, which is not included here, shows that this equation can be 
replaced by the locally equivalent equation 



l + logp W (x,p)+0((logp) p ) =0, 



(3.13) 
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where p is a positive number, p > 1. Then, from Lemma 

1 ..... .. .,, ^ 3c*(l,l) .^.^^ _,,, Tr ,_ ,,„„.. ,,, 



log(^)c*(l,l) ^+log(px)$(x,p) + *(x,p) + 0( a; - d ), (3.14) 

which implies p. lip . 

For the computation of a'(l,p), one uses the simple relation 

pa (l,p) = ai, 

a 

obtained by derivation of (|2.2|) at z = 0, and (|3.12p follows. □ 

4 Remarks and complements 

It is important to note that one of the main technical difficulties of the problem, 
besides its strongly non-linear feature, comes from the phase transition which 
appears for p = 1. Actually, c(z,p) > for < z < and c(£(p),p) = 0. 
Then, for z 3> c(z, p) has wild unbounded oscillations. In particular, this 
implies that L c(x, p)dx does not exist. 

On the other hand, when p = 1, the c(z, 1) is no more an integral function: 
it has a singularity (which seems to be a pole of order 3) located on the negative 
real axis; it does not vanish for z > 0, and the quantity J Q c(x, l)dx = A is 
finite. In the half plane 5R(z) > 0, the following expansions hold: 

c*(l,l), 9 , Dloez /log 2; 

c(z,l) = exp )-J—Llog 2 z + Blogz+ ^- + o 1 

2 z z 

v(z,l) = c^l^logz + B + ^+o^ 

z l 

where B and D are some constants. 

Finally, iterating (|3.2p . one could improve some of the estimates given in the 
previous section, rewriting c(z, p) as 

c(z,p) = c(zp I+1 ,p)Y[[l + log pv(zp\p)], 

i=0 

where / is an arbitrary positive integer. The above product is uniformly conver- 
gent, V/ < 00, for all z in a compact set of the complex plane, since it behaves 
like the series 

z(logp) 3 £5>*+ fe +M^ +fc+ \p). 

fe>0 i>l 

This series has its modulus bounded by \zp 2 c{pz, p)\, so that, from the max- 
imum modulus principle, it converges uniformly. Then, one can choose I to 
ensure 

, log z 

< zp 1 < 1, i.e. I < -—2-. 

logp 

and make use of the properties of the T function to estimate the product. 
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